Ground-based multisite observations of two transits of HD 80606b 
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ABSTRACT 

We present ground-based optical observations of the September 2009 and January 2010 tran- 
sits of HD 80606b. Based on 3 partial light curves of the September 2009 event, we derive a 
midtransit time of T c [HJD] = 2455099.196 ± 0.026, which is about la away from the previously 
predicted time. We observed the January 2010 event from 9 different locations, with most phases 
of the transit being observed by at least 3 different teams. We determine a midtransit time 
of T c [HJD] = 2455210.6502 ± 0.0064, which is within 1.3er of the time derived from a Spitzer 
observation of the same event. 

Subject headings: planetary systems — stars: individual (HD 80606) 
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1. Introduction 

The sample of transiting exoplanets has grown 
rapidly in recent years, but HD 80606b stands 
out from this crowd by virtue of its long pe- 
riod (111.4 days) and high orbital eccentricity 
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(e=0.93). Initially the planet was discovered in 
a Doppler survey and was not known to tran- 
sit (Naef et al. 2001). However, the planet is 
near pericenter during superior conjunctions, in- 
creasing the probability of occultations (secondary 
eclipses) . T his motivated the observations of 



Laughlin et al.l (j2009l ). who found that the or- 
bital inclination is indeed close enough to 90° for 
occultations to occur. Soon after, the discovery 
of a transit (p rimary eclipse) was reported by 
several authors dMoutou et al.l 2009; iFossev et al 



20091 : iGarcia-Melendo fc McCulloughl 12009^ In 



addition, the orbit was shown to be misaligned 
with the plane defined by the stellar rotation 



(Moutou et al. 2009; 


Pont et al. 


20091 fWinn et al. 


20091; Hebrard et al. 


2010). This system has be- 



study for theories of 
and giant 



come an important case 
orbital migration, tidal in teractions 
planet atmospheres ( e.g..^Wu fc Murr m 



Fabrvckv fc Tremainel 120071 ; Triaud et al 



Pont et al. 2008; La ughlin et al.l2009HKnutson et al 



2003; 



2010; 



20091) . 

For any transiting planet, it is desirable to ob- 
serve multiple transits. This allows the system 
parameters to be refined, especially the orbital pe- 
riod. In addition, a deviation from strict period- 
icity may be the signature of additional bodies in 
the s ystem, such as another p l anet or a satellite 
(e.g.. iHolman fc Murravl l2005t lAeol et al.l l2005t 
ISimon et alJl2007t iNesvornv fc Beauge1l201oj l~ 

In the case of HD 80606b observations of com- 
plete transits are challenging. This is partly be- 
cause transits are rare, occurring once every 111.4 
days, and also because the transit duration is 
nearly 12 hours, making it impossible to be ob- 
served in en t irety from a single ground-ba sed site. 
Hidas et al.l ((2010) and I Winn etail (|2009l) carried 



out multisite campaigns to achieve more complete 
coverage of the transit. This paper reports the re- 
sults of our more recent observations. In Septem- 
ber 2009, only the first portion of the transit could 
be observed, but our combined light curve for Jan- 
uary 2010 ranges over the entire transit, with most 
phases having been observed with at least 3 dif- 
ferent telescopes. 

Hebrard et al. I (|2010h utilized the Spitzer obser- 



vatory to obtain a high quality light curve of the 
entire January 2010 transit event. Besides deriv- 
ing refined system parameters, they measured a 
midtransit time earlier by 3cr than was predicted 



by I Winn et al. ( 2009h . This discrepancy suggested 
the intriguing possibility of the existence of a third 
body in the HD 80606 transiting system. The pri- 
mary goal of our analysis was to determine the 
times of these events as accurately as possible. Al- 
though our results are not as precise as the Spitzer 
result, we provide an independent estimate of the 
midtransit time. Our observations and photomet- 
ric processing are described in § [2j In § [3] we de- 
scribe our analysis, the results of which are pre- 
sented in § |4] and discussed in § [5] 

2. Observations 

The new data presented in this paper are 3 
photometric time series spanning the first portion 
of the transit of 2009 September 23/24 (hereafter 
Sep09), based on observations across North Amer- 
ica; and 10 photometric time series spanning the 
entire transit of 2010 January 13/14 (hereafter 
JanlO), based on data from Israel, Europe, Ca- 
nary Islands, North America and Hawaii. 

In all cases we gathered CCD images encom- 
passing the target star, HD 80606 (G5V, 7=9.06, 
B — 7=0.76), and its visual binary companion, 
HD 80607 (G5V, 7=9.07, B - 7=0.87) which 
is located 21" east. The similarity in brightness 
and color between the two stars facilitates differ- 
ential photometry. Since the stars are relatively 
bright, we defocused most of the telescopes that 
were used, thereby allowing for longer exposures 
without saturation, and reducing the impact of 
pixel sensitivity variations and seeing variations. 
We always ensured that the PSFs of the two stars 
were well separated. Brief descriptions of the ob- 
servations from each site are given below, along 
with an abbreviated observatory name that we will 
use in the remainder of this paper. 

Wise Observator$\ (WO), Israel. The target 
was observed with the 1.0 m telescope for 9 hours 
on the JanlO transit night, completely covering 
ingress, and for 2 hours on the following night. 
An RGO Z filter was used and the telescope 
was slightly defocused. The detector was a back- 
illuminated Princeton Instruments CCD, with a 
12/6 x 13/0 FOV and a pixel scale of 0."580 pixel" 1 . 

Gran Telescopio Canaria^ (GTC), La Palma, 



1 http://wise-obs.tau.ac.il/ 

2 http : //www .gtc . iac . es/en/pages/gtc .php 
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Canary Islands. We observed the target with the 
Optical System for Imaging and low Resolution In- 
tegrated Spectroscopy (OSIRIS), mounted on the 
10.4 m GTC telescope. The observations lasted for 
8.3 hours on the JanlO transit night, all within the 
transit. A small defocus was applied, resulting in 
a typical PSF FWHM of 1.5". We used the narrow 
band imager with the red range tunable filter, al- 
ternating between four narrow bands in the range 
7680-7780 A, each with a FWHM of 12 A. The in- 
strument has a maximum FOV of 7/ 8 x 7/ 8 and a 
pixel scale of 0."127 pixel -1 . Analysis of the light 
curve for ea ch individua l narr ow band filter are 
presented in IColon et all (|2010h : here, we are con- 
cerned with the light curve based on a summation 
of the flux observed in all four filters. We adopt 
the SDSS-i' band limb darkening coefficients when 
fitting this light curve. 

Observatoire de Haute Provenc^ ( OHP), France. 
The target was observed with the 1.2 m tele- 
scope for 4.3 hours during the first half of the 
JanlO transit, with no out-of-transit observations. 
As was the case with the spectroscopic observa- 
tions secured simulta neously with Sophie at OHP 
( Hebrard et al. 2010l ). the transit sequence had to 
be stopped due to cloudy weather. Photometric 
observations were also secured on several nights 
shortly before and after the transit night, to put 
constraints on the stellar a ctivity. These data ar e 
described in more details bv lHebrard et al.l (|2010l ). 
No defocus was applied here; instead, to allow for 
longer exposures, we used a neutral density filter 
along with a Gunn-r filter. 

Allegheny Observator^ (AO), Pittsburgh, Penn- 
sylvania, USA. We used the 0.41 m (16 in) Meade 
telescope and a SBIG 0."56 pixel -1 CCD, with 
a FOV of 20/0 x 30/0. Observations of the 
JanlO transit were done in the / filter with no 
defocus, from the time of second contact until 
about 2 hours after the transit ended, a total of 
11.4 hours. 

Rosemary Hill Observator^ (RHO), Bron- 
son, Florida, USA. We used the 0.76 m Tinsley 
telescope and a SBIG ST-402ME CCD camera 
mounted at the f/4 Newtonian focus, with a FOV 
of 3/88 x 2/56 (half of the entire FOV). During 



the JanlO event we observed for 4.9 hours with no 
out-of-transit data. We used the same instrument 
to observe the Sep09 event, obtaining 2.4 hours of 
data just before the transit started, which helps 
constrain the transit start time. Both events were 
observed with the SDSS-i' filter and the telescope 
was defocused. 

Fred Lawrence Whipple Observator^ (FLWO), 
Mt. Hopkins, Arizona, USA. We used Keplercam, 
mounted on the 1.2 m telescope, and observed the 
target for 2.8 hours during the flat bottom part of 
the JanlO transit. The same instrument was also 
used to observe the very beginning of the Sep09 
event, obtaining 1.8 hours of data. KeplerCam 
includes a 4K x 4K FairChild Imaging CCD486, 
with a pixel scale of 0."672 pixel -1 (2x2 bin- 
ning) and a 23/1 x 23/1 FOV. FLWO observations 
of both transit events used the SDSS-i' filter and 
started immediately after the target rose. 

Table Mountain Observator$\ (TMO), Wright- 
wood, California, USA. We used the 1.0 m tele- 
scope with an Apogee U16 4K x 4K CCD, with 
a pixel scale of 0."18 pixel -1 and a 12/3 x 12/3 
FOV. The SDSS-i' filter was used here. TMO ob- 
servations completely cover the JanlO egress, with 
a total of 6.6 hours, of which 2.9 hours are after 
fourth contact. 

George R. Wallace Jr. Astrophysical Observa- 
tor$ (WAO), Westford, Massachusetts, USA. We 
used two identical 0.36 m telescopes on the JanlO 
transit night, each equipped with a SBIG STL- 
1001E IK x IK CCD camera with a pixel scale of 
l."29 pixel -1 , resulting in a FOV of 21/5 x 21/5. 
We refer to these instruments as WAO-1 and 
WAO-2. At both telescopes we observed in the I 
filter with no defocus. With WAO-1 we were able 
to observe for 15 minutes during the flat bottom 
part, before observing conditions degraded. Ob- 
servations resumed with both instruments at the 
time of third contact and continued for 4.2 hours, 
of which 2.4 hours were out-of-transit. 

Faulkes Telescope Nortf^ (FTN), Mt. Haleakala, 
Maui, Hawaii, USA. We used the LCOGT0 2.0 m 
FTN telescope, and the Spectral Instruments cam- 



3 http://www.obs-hp.fr 

4 http : //www.pitt . edu/$\sim$aobsvtry/ 

5 http : //www . astro . uf 1 . edu/ inf ormat ion/rho . html 



6 http : //www . sao . ar izona . edu/FLWQ/whipple . html 

7 http://tmoa.jpl.nasa.gov/ 

8 http://web.mit.edu/wallace/ 

9 http://faulkes-telescope.com/ 
10 http://lcogt.net 
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era with the Pan-STARRS Z filter. The camera 
has a back-illuminated Fairchild Imaging CCD 
and we used the default 2x2 pixel binning mode, 
with an effective pixel scale of 0."304 pixel . The 
telescope was defocused and the 10.'5 x 10.'5 FOV 
was positioned and rotated so the guiding cam- 
era FOV will contain a suitable guide star. We 
observed the target on the JanlO transit night 
for 9.2 hours, from the beginning of egress until 
7 hours after the transit ended. We observed also 
during the two adjacent nights, for 2 hours on the 
preceding night (Jan. 12/13), and for 3.5 hours on 
the following night (Jan. 14/15). 

Mount Laguna Observator^\ (MLO), San 
Diego, CA. Only the Sep09 event was observed. 
We used the 1.0 m telescope with a Fairchild Imag- 
ing 2K x 2K CCD and SDSS-i' filter. Observa- 
tions were done with a 300 x 300 pixel sub-array 
and a FOV of 2/0 x 2/0. The telescope was de- 
focused and the very beginning of the event was 
observed, for 1.6 hours with no out-of-transit data. 

Four additional light curves were obtained 
for the JanlO event by the Liverpool Telescope, 
MONET-North telescope, DeKalb Observatory, 
and a third telescope at WAO, and one for the 
Sep09 event by AO. However, those data displayed 
a very high noise level and strong systematic ef- 
fects, and were not included in our subsequent 
analysis. 

The CCD data were reduced using standard 
routines for bias subtraction, dark current sub- 
traction (when necessary) and flat-field correction. 
We used aperture photometry to derive the flux of 
HD 80606 and HD 80607, and divided the former 
by the latter to obtain a time history of the flux 
ratio, which we refer to as the light curve. We 
took care to choose aperture sizes to avoid con- 
tamination of one stellar signal by the other star. 
Our time stamps represent the Heliocentric Julian 
Date, based on the UTC at midexposure (and not 
the unifo rmly flowing ter r estria l time system advo- 



cated by lEastman et ah (2010)). This is also the 



ti me system that was u sed for the Spitzer analysis 



of Hcbrard et al 



(|20inh . 

All light curves were averaged into 10 minute 
bins, using 3tr outlier rejection. The error bar as- 
signed to each data point was the standard devi- 
ation of the mean of all the measurements con- 



11 http://mintaka.sdsu.edu/ 



tributing to each 10 minute bin, which ranged in 
number from 4 to 63. There is no significant in- 
formation loss due to time binning, because the 
binning time of 10 minutes is shorter than the du- 
ration of ingress and egress by more than an order 
of magnitude. 

The data cover only the first portion of the 
2009 September 23/24 event, but they provide 
complete coverage of the 2010 January 13/14 
event. Although the first hour of the transit was 
observed from only one observatory (WO), the rest 
of the transit was observed by 3-5 different sites, 
which is very helpful for identifying and decreas- 
ing the influence of any systematic effects that are 
specific to each observatory (i.e. correlated noise, 
or red noise) which is frequently a pr oblem with 
ground-based photometric data (e.g.. IPont et al.l 
l2006t ICarter fc Winnll2009HSvbilski et al.ll2010h 

Table Q] gives the photometric data that were 
obtained and analyzed. Each data point repre- 
sents a 10 minute binned average of the flux ratio 
of HD 80606 to HD 80607. The normalization fac- 
tors and error rescaling factors that are described 
in § [3] were not applied to the data given in the 
table. Table [5] gives a list of all the observatories. 

3. Data Analysis 

In this section we describe the methods by 
which we combined the data and derived the mid- 
transit time of each event. We describe the process 
in detail for the JanlO event; the details were very 
similar for the Sep09 data. 

Because the quality of the Spitzer light curve 
is superior to any ground-based light curve, we 
did not attempt to use our data to refine the ba- 
sic system parameters other than the midtransit 
ti mes. Instead, we use d the parameters derived 
bv lHebrard et all (|2010f ) as constraints on the light 
curve shape, while allowing the midtransit time to 
be a free parameter, as described below. 

A simultaneous analysis of several light curves 
obtai ned by different ins trume nts is a chal l enging 
task. IWinn et al.l (|2009l ) and Iffidas et al.l (|2010l) 
have carried out a similar task, although for a 
smaller number of data sets. One of the crucial 
points is the placement of the different flux ratio 
light curves onto the same scale, despite the differ- 
ences in bandpasses, detectors, and weather con- 
ditions at each observatory. One way of thinking 
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about this problem is that we need to establish 
the out-of-transit flux ratio that was measured, 
or that wou l d be measured, by each observatory. 
Winn et ah! (2009) performed this calibration by 
using data taken on nights when the planet was 
not transiting, a method that may be affected by 
night-to-nig ht variations due to varying observing 
conditions. Hidas et al. ( 2010h allowed the over- 
all flux ratio scale to be a free parameter for each 
light curve, thereby increasing the overall number 
of fitted parameters. Here we chose an intermedi- 
ate approach. We assigned a normalization factor 
to each light curve, estimating it from the out-of- 
transit data whenever possible, and allowing it to 
be a free parameter when there was insufficient 
out-of-transit data. 

We divided our light curves into three groups, 
based on the amount of out-of-transit data: 

Group I includes the light curves with abun- 
dant out-of-transit information. The only mem- 
bers of this group are the WO and FTN light 
curves of the JanlO event, where out-of-transit 
measurements were obtained on the transit night 
and on one (for WO) or two (for FTN) of the ad- 
jacent nights. For each of the two light curves 
separately, we subsequently fitted a 2nd degree 
polynomial to the out-of-transit flux ratio points 
vs. airmass, time and PSF FWHM, and divided 
the entire light curve by this polynomial. The ef- 
fect of this process on the in-transit points was 
small, typically at the few 0.01% level. The as- 
signed normalization factor for each of the two re- 
sulting light curves was 1.0. 

Group II includes light curves that have at least 
1.5 hours of out-of-transit measurements. Their 
normalization factors were taken to be the mean 
out-of-transit flux ratio. This group includes 4 
light curves: the AO, TMO, and the two WAO 
light curves. 

Group III includes the remaining 4 light curves 
(GTC, OHP, RHO and FLWO) with either a small 
amount of data or no out-of-transit data at all. 
The normalization factors were taken to be free 
parameters in our model. 

Our model for the data is based on the premise 
of two spherical objects, a non-luminous planet 
and a limb-darkened star with a quadratic limb 
darkening law, in an eccentric Keplerian orbit. 
For each binned time stamp we calculated the 



sky projecte d planet-star dis t ance and used the 
equations of iMandel k, Ago! (|2002l ) to calculate 
the relative flux at that time. Our code ac- 
co unts for the light tra vel time effects described 
by IHebrard et al. I (l2010h . 

The model included a total of 34 parameters: 
orbital period P, planet-to-star radius ratio r = 
R p l R S1 orbital scmimajor axis in units of the stel- 
lar radius a/R s , orbital eccentricity e, argument of 
periastron to (in fact our fitting parameters were 
actually ecosw and esinw), inclination angle i, 
an individual periastron passage time {Tp^}^ for 
each of the 10 light curves (which were later con- 
verted into midtransit times), two limb darkening 
coefficients u\ and U2, for each of the different 4 
filters used, and 10 normalization factors, one for 
each light curve. 

The limb darke ning co e fficien t s wer e estimated 
from the grids of IClaretl (|2000l |2004| ) for a star 
with T fiff =5645 K, log£=4.5 and [Fe/H]=0.43 
(jNaef et al.l 120011 ). and were held fixed in the fit- 
ting process as the data are insensitive to the 
coefficients. As mentioned earlier, we used the 
parameters of IHebrard et al.l ( 20101 ) to constrain 
the light curve shape. We used their values and 
uncertainties for P, r, a/R Sl ecosw, esinw and i 
as a priori Gaussian constraints by adding penalty 
terms to the % 2 fitting statistics. We constrained 
in a similar way the normalization factors of the 
6 light curves in Groups I and II described above, 
while assuming a normalization factor uncertainty 
of 0.1%. This uncertainty is larger than that of 
the light curves mean out-of-transit flux ratio, 
typically a few times 0.01%, and it was used in 
order to prevent the fitting process from being 
dominated by a single light curve. 

Out of the 34 model parameters, 8 were held 
fixed (the limb darkening coefficients), 12 were 
controlled mainly by Gaussian priors (P, r, a/R s , 
ecosw, esinw, i and 6 normalization factors), and 
the remaining 14 were free parameters with uni- 
form priors ({T Pi i}}^i and 4 normalization fac- 
tors). Our fitting statistic was: 



X 



2,2 I 2 



(1) 



where the first term on the right hand side is the 
usual x 2 statistics: 



431 



4=E 



i=l 



fi(obs) - fi(model) 



(2) 
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the second term includes the penalties for the light 
curv e parameters (values and uncertainties taken 
from iHebrard et all I2010L P in days and i in de- 
grees): 



Xorb 



P - 111.4367 



0.0004 



a/R s - 97.0 



n 2 



1.6 



-0.1001 
0.0006 

ecosu - 0.4774 
0.0018 



esinw- (-0.8016)" 


2 

+ 


~i- 89.269" 


0.0017 




0.018 



+ 
+ 

2 



(3) 



and the third term constrained the normalization 
factors of the 6 light curves for which we have suf- 
ficient out-of-transit data: 



X 2 

A-norrn 



E 

1=1 



f, 



oot,l 



oot.l 



0.001 



(4) 



After a preliminary fit, the residuals of each 
data set were carefully examined. Two data points 
were clear outliers, departing from the model by 
>4cr, and were rejected, leaving a total of 431 
points. The rejected points were either the first 
or last data points in the time series, and were 
probably affected by the high airmass or relatively 
bright twilight sky. After refitting there were no 
additional >4er outliers. 

Next, for the purpose of determining parameter 
uncertainties, we determined appropriate weights 
for each data point. This was done in two steps 
specific to each light curve. First, we rescaled 
the error bars such that the median error bar 
was equal to the rms residual. We named this 
rescaling factor a. If the median error bar was al- 
ready equal to or smaller than the rms residual, 
we set a = 1.0. Second, we attempted to account 
for correlated (red) noise in each time series on 
the critical time scale of the ingress/egress dura- 
tion (2.8 hour s). We used the "time-averaging" 
method (e.g., Winn et al. 20081) . in which the 



residuals are binned using several bin sizes, close 
to the duration of ingress and egress. The amount 
of correlated noise is then quantified by the ratio 
between the binned residual light curves standard 
deviation and the expected standard deviation as- 
suming pure white noise. For each light curve we 
took j3 to be the largest ratio among the bin sizes 
used, and multiplied the individual error bars by 



that factor. We took /3 to be 1.0 when the time- 
averaging method gave a smaller value. The val- 
ues of a and j3 for each light curve are listed on 
Tabled 

To determine the "best" values of the param- 
eters, and their uncertainties, we used a Monte 
Carlo Markov Chain (MCMC) algorithm (e.g., 
iTegmark et al.ll2004l : lFordll2005h with Metropolis- 
Hastings sampling, which has become the stan- 
dard practice in the literat ure on transit photom- 
etry " ' ' ' V /■<<:./■. . , 1 

2007 



(Hol man et all 2006 
Burke et alJl2007t ) 



Collier Cameron et al 



We used here an adap- 
tive approach , simil ar to the one described by 
Shporer et all (|2009h . The algorithm steps from 



a multidimensional point in parameter space, Pi, 
to the next, Pi+i, according to 



Pi+i = Pi + f*G(0, 1), 



(5) 



where G(0, 1) is a vector of numbers picked ran- 
domly from a Gaussian distribution of zero mean 
and unit standard deviation, a is a vector of the 
so-called step sizes, and / (the only scalar in 
Eq. [5]) is a factor chosen to control the fraction 
of accepted steps. The value of / was readjusted 
every 10 3 steps, to keep the acceptance fraction 
near 25% (|Gregory|[2005h . Our final MCMC in- 
cluded 10 long chains of 500,000 steps each, start- 
ing from different initial positions in parameter 
space spaced apart by ~5cr from the best-fitting 
parameters. The posterior probability distribu- 
tion of each parameter was constructed from all 
long chains after discarding the first 20% of the 
steps. We took the distribution median to be the 
"best" value and the values at the 84.13 and 15.87 
percentiles to be the +la and -la confidence un- 
certainties, respectively. 

4. Results 

Results of the fitting process are presented vi- 
sually in Fig. Q] where each light curve is plotted 
separately and overplotted by the fitted model, us- 
ing limb darkening coefficients of the correspond- 
ing filter. Fig. [5] shows the combined light curve, 
with the 13 JanlO and Sep09 light curves plotted 
on top of each other. 

The fitted model for the JanlO event has 



442/429, and for the Sep09 event 



X /Ndoi = 36/33. The rms residual of each binned 
light curve is given in the second-to-last column 
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of Table H The rightmost column of Tabled 
lists the photometric noise rate (PNR) of the un- 
binned light curve, defined as the rms divided 
by VT, where T is the median number of data 
points per unit time (the "cadence"). The PNR 
is meant to be a quantitative comparison between 
the statistical pow er o f different data sets (se e also 
Burke et al.ll2008l and IShporer et al~ll20Qc Ph . 



The fitting process resulted in 10 estimates of 
the periastron passage time, {T Pj j}J ; £ 1 . Using 
the resulting distributions of the light curve pa- 
rameters we numerically translated each perias- 
tron time to midtransit time. We then averaged 
the midtransit times to get our final estimate of 
the JanlO midtransit time, T Cl j an io. While ex- 
amining the 10 individual midtransit time esti- 
mates we noticed that 3 of them have relatively 
large uncertainties, larger than 0.015 days, while 
the rest have uncertainties in the range of 0.005- 
0.010 days. Therefore we removed them from the 
sample before averaging. Those 3 are the OHP, 
RHO and FLWO light curves, spanning less than 
5 hours without any out-of-transit data, hence the 
increased uncertainties are expected. We used 
an unweighted average to get our final midtransit 
time estimate, of T CiJanW [HJD] = 2455210.6502, 
with an rms of 0.0064 days, close to the typical 
uncertainty of the 7 individual estimates. Using 
regular averaging acts to average out possible cor- 
related noise affecting individual estimates. Using 
median or weighted average changes the midtran- 
sit time by < 0.25<r. Including the 3 large un- 
certainty T c /s mentioned above results also in a 
small change to T Ci j an xo, of less than 0.3cr but the 
scatter is increased by 50%. 

To check the sensitivity of our result to our 
treatment of correlated noise, we reran our anal- 
ysis for the JanlO data while using /3 = 1.0 for 
all light curves. The resulting T c ^ uncertainties 
were smaller by typically 25%-30%, and the aver- 
age midtransit time was 2455210.6476 ± 0.0093, 
earlier by 0.4er than our preferred analysis. 

For estimating the midtransit time of the Sep09 
event, T Ci s e p09, we used only 2 of the 3 T C) j's 
as the RHO data is entirely out-of-transit. The 
average between these two is T Ct s ep o9[HJD] = 



1 2 We note tha t the equation for the PNR given by 
IShporer et all | J200gf l includes a typographical error, al- 
though their calculations are correct 



2455099.196 and we assign it an uncertainty of 
0.026 days, similar to that of the two T Cl j's. 

We note that figs. [Hand [2] were produced using 
the average midtransit times of each event, derived 
above, not with the T Cj j of each partial light curve. 

5. Discussion 

Table [3] lists the mi dtransit times derived here, 
the one measured by Hebrard et al. ( 2010l) and 
those predict ed by I Winn et al.l (120091) . For the 
JanlO event Hebrard et al. ( 2010l ) measured a 
Tr...TanAo wh i ch is 3(7 earlier than predicted by 
Winn et al. (2009). Our m easured T r ...j nn m is in- 



ter mediate betw e en tho se of Hebrard et al. ( 201dh 



and Winn et al 



Hebrard et al.l 



(200$. It is 1.3cr later than the 



2010) time, and earlier by 1.1a 



than the Winn et al. (2009) prediction. There- 
fore, it does not confir m nor refute the ear lier 
transit time measured by IHebrard et al. ( 2010h . 



Th ere is s ome discrepancy be tween lHebrard et al 
(|2010l ) and IWinn et al.l (|2009h also in the three 
parameters determining the light curve shape, 
namely r, a/R s and i. Specifically, the values 
for the planet-star radius ratio differ by about 
2.5cr, and for a/R s and i the difference is at the 
1.5f7 level. To check which set of values is pre- 
ferred by ou r data we reran our analysis using the 
Winn et al.1 (|2009h values for these three param- 
eters as prior constraints in Eq. [3] This resulted 
in an increased Y 2 /-/Vdof of 40% showi ng that 
our data prefers the IHebrard et al. ( 2010l ) values. 
This indicates the discrepancy is not the result 
of a wavelength dependent radius ratio and more 
likely results from underestimated uncertainties. 

IHebrard et al.l (|2010l) noticed a "bump ' in their 
Spitzer light curve: a ~0.1% increase in flux last- 
ing about 1 hour, just before midtransit. This 
bump could be caused by the passage of the 
planet in front of a dark spot on the star, al- 
thou gh HP 80606 is not known to be an active 
star (He brard et al. 2010l) . The decreased sur- 



face brightness of a spotted surface element re- 
sults from decreased temperature relative to a 
non-spotted surface element, so the flux increase 
during a spot crossing is expected to be wave- 
length dependent, increasing with shorter wave- 
lengths (|Pont et alj|2008t iRabus et al.ll2009T) . 

The inset in Fig. [2] shows that our data are 
not sensitive enough to identify a 0.1% flux vari- 
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ation, although it does show that a 0.2% in- 
crease is unlikely. We note that another phe- 
nomenon that could, in principle, result in a 
small brief wavelength-independent bump is a 
triple conjunction of the s tar, planet and a moon 



increase in t ransit timing variations with orbital 



orbit i ng the planet (e.g ., ISartoretti fc Schneider 



19991 : Simon et al. 20101 ). A moon orbiting the 



planet could also be responsible for a shifted mid- 



transit 



Simon et al 



time (e.g ., 
2007). 



Sartoretti fc Schneider 19991 : 



Although the existence of a 
moon is an exciting possibility, we caution that it 
is not likely to be the case in the current system 
due to the improbability of a triple conjunction, 
and perhaps also the dynamical instability of such 
a scenario. 

We have presented here the results of an obser- 
vational multisite campaign, along with a method 
to combine the partial ground-based light curves 
to obtain a complete transit light curve of a planet 
with a long period and long transit duration. 
Our results in comparis on with the simultan eous 
Spitzer observations of Hebrard et al. ( 2010() al- 
lows an assessment of the quality of multisite cam- 
paigns and develop this method for the future dis- 
coveries of additional systems with similar charac- 
teristics. Such systems will be discovered by pho- 
tometric observations of stars known to have plan- 
ets from RV surveys, like HD 80606 itself, and by 
the spaced-based transiting planet hunters C0R0T 
and Kepler, capable of continuous photometric 
monitoring. The recent discovery of CoRoT-9b 
(|Deeg et al.ll2010l) . with a 95 day period and an 8 
hour long transit, is an excellent example. Only 
two transit events were observed by C0R0T, al- 
lowing accurate determination of the light curve 
shape parameters, while ground-based photomet- 
ric follow-up observations were needed to refine 
the transit ephemeris. Using accurate light curve 
parameters, determined by space-based observa- 
tions, to measure a midtransit time of a ground- 
based light curve is similar to the approach we 
have taken here. 

Ground-based observations are limited by ef- 
fects of the atmosphere and short observing win- 
dows (the latter leading to the need to combine 
data obtained by different instruments) , but small 
to medium telescopes are relatively easily accessi- 
ble, compared to space telescopes. Although form- 
ing a collaboration between several observatories is 
not a trivial task, those will be motivated by the 



perio d (e.g., iHolman fc Murravl l2005t lAgol et al 
120051 ). thus allowing ground-based observations to 
detect the effects of an additional, unseen planet 
in the system. 

Some of the limitations and difficulties men- 
tioned above can be minimized by using a network 
of identical instruments spread around the globe, 
allowing for continuous (24/7) monitoring of the 
target, especially when the observing windows 
partially overlap, allowing an accurate combina- 
tion of the partial light curves . The LCOGT net - 
work (e.g., Brown et al. 2010t Lewis et al. 2010h . 
once completed, is meant to be such a network, 
in both the Northern and Southern hemispheres. 
For targets in Northern (or Southern) positions 
such as HD 80606, the observational windows will 
partially overlap. 
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Table 1 

hd 80606 light curves presented in this work 1 



HJD Flux Ratio Error Observatory 2 



January 2010 transit event 



2455210.288327 


1.00114 


0.00116 


WO 


2455210.445392 


1.11691 


0.00024 


GTC 


2455210.446272 


0.99860 


0.00106 


OHP 


2455210.515999 


1.10749 


0.00406 


AO 


2455210.607419 


1.11228 


0.00107 


RHO 


2455210.654098 


1.11259 


0.00110 


FLWO 


2455210.748739 


1.11077 


0.00142 


TMO 


2455210.588711 


1.10630 


0.00116 


WAO-1 


2455210.815045 


1.11449 


0.00134 


WAO-2 


2455209.863305 


0.99977 


0.00037 


FTN 


September 2009 transit event 


2455098.867697 


1.12197 


0.00227 


RHO 


2455098.950735 


1.11411 


0.00138 


FLWO 


2455098.974510 


1.12694 


0.00101 


MLO 



1 The table is given in its entirety in the on-line version 
of the manuscript. A sample is given here, to show its 
format. 

2 Observatories name abbreviations: WO = Wise Obser- 
vatory, GTC = Grand Telescopio Canarias, OHP = Ob- 
servatoire de Haute Provence, AO = Allegheny Observa- 
tory, RHO = Rosemary Hill Observatory, FLWO = Fred 
Lawrence Whipple Observatory, TMO = Table Mountain 
Observatory, WAO = Wallace Astrophysical Observatory, 
FTN = Faulkes Telescope North, MLO = Mountain La- 
guna Observatory 
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Table 2 
Observatories List 



# 


Obs. 


Tel. 


Filter Transit Part 1 




a 


13 


RMS 2 


PNR 3 






[m] 


OIBEO 








[%] 


[% min- 1 ] 




Jan 2010 


1 


WO 


1.0 


Z O I B - 




1.45 


1.00 


0.11 


0.29 


2 


GTC 


10.4 


7680-7780A (j!) - I B E 




3.29 


2.15 


0.05 


0.09 


3 


OHP 


1.2 


Gunn-r+ND 4 - I B - 




2.49 


2.27 


0.13 


0.25 


4 


AO 


0.41 


/ --BE 





1.00 


1.10 


0.12 


0.43 


5 


RHO 


0.76 


i' --BE 




1.91 


2.59 


0.09 


0.15 


6 


FLWO 


1.2 


i' - - B - 




1.11 


1.53 


0.08 


0.21 


7 


TMO 


1.0 


i' --BE 





1.04 


1.18 


0.07 


0.20 


8 


WAO-1 


0.36 5 


I --BE 





1.64 


1.00 


0.16 


0.33 


9 


WAO-2 


0.36 5 


/ E 





1.25 


1.35 


0.11 


0.28 


10 


FTN 


2.0 


Z E 





1.24 


1.07 


0.07 


0.23 


Sep 2009 


1 


RHO 


0.76 


i' - - - 




1.00 


1.00 


0.12 


0.44 


2 


FLWO 


1.2 


i' o I - - 




1.00 


1.00 


0.08 


0.30 


3 


MLO 


1.0 


i' o I - - 




1.06 


1.00 


0.09 


0.23 



1 OIBEO for Out of transit before ingress, Ingress, flat Bottom, Egress and Out of transit after egress. 
2 rms residual of the 10 minute binned light curve. 

3 Photometric Noise Rate of the unbinned light curve, calculated as rms/VT where T is the median number 
of exposures per minute. 

4 Neutral Density filter. 

5 Two identical telescopes. 



11 



Table 3 

MlDTRANSIT TIMES. 



Reference T c [HJD] 

Jan 2010 

This work 2455210.6502(64) 
Hebrard et al. 2010 2455210.6420(10) 
Winn et al. 2009 1 2455210.6590(51) 

Sep 2009 

This work 2455099.196(26) 
Winn et al. 2009 1 2455099.2216(50) 



1 Predicted T c based on the ephemeris 
given in that reference. 



Winn, J. N., et al. 2009, ApJ, 703, 2091 
Wu, Y., & Murray, N. 2003, ApJ, 589, 605 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 
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10.1 10.2 10.3 10.4 10.5 10.6 10.7 10.8 10.9 11 11.1 

HJD - 2455200 



1.01 - 




98.7 98.8 98.9 99 99.1 99.2 99.3 99.4 99.5 99.6 99.7 

HJD - 2455000 



Fig. 1. — New light curves of HD 80606. Top. — The 10 light curves of the JanlO event. Bottom. — The 3 light 
curves of the Sep09 event. Each light curve is overplotted by the best-fitting model, with limb-darkening 
coefficients appropriate for the observing bandpass. 
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1 .005 - 




-10 -5 5 10 

Time from mid transit [hours] 



Fig. 2. — Composite light curve. The phase-folded light curve based on the 13 new data sets analyzed in 
this paper. Light curves taken with different filters are represented by different colors. For each filter the 
corresponding model light curve is plotted as a solid line with the same color. Data obtained in the Z band 
is in gray, i and / data in black and r data in red. Open circles represent the Sep09 data and filled circles 
the JanlO data. Plotted below the light curve are the residuals, with error bars. The dashed line marks 
the residual zero point, and the 3 dotted lines mark relative flux residual levels of 0.1%, 0.2% and 0.3%. 
The inset shows a zoome d-in view of th e resid uals during the phase when a "rebrightening" or "bump" was 
observed with Spitzer bv lHebrard et all (|2010h . The Spitzer feature is overplotted in blue. Our data neither 
confirm nor refute the existence of the bump. 
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